Towards a Landau-Zener formula for an interacting Bose-Einstein condensate 
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We consider the Landau-Zener problem for a Bose-Einstein condensate in a linearly varying two- 
level system, for the full many-particle system as well and in the mean-field approximation. The 
many-particle problem can be solved approximately within an independent crossings approximation, 
which yields an explicit Landau-Zener formula. 
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I. INTRODUCTION 

During the last years, a lot of work has been de- 
voted to the nonhnear Landau-Zener problem, which 
describes a Bose-Einstein condensate (BEG) in a time- 
dependent two-state system in the mean-field approxima- 
tion 0,0, 13- As in the celebrated original Landau-Zener 
scenario, the energy difference between the two levels is 
assumed to vary linearly in time. This situation arises, 
e.g., for a BEG in a double- well trap or for a BEG in an 
accelerated lattice around the edge of the Brillouin zone. 
A major question in such a situation is the following: Ini- 
tially the two states are energetically well separated and 
the total population is in the lower state. Then the en- 
ergy difference varies linearly in time, such that the two 
levels (anti-) cross. Finally the states are energetically 
well separated again, however they are just exchanged. 
What is the probability of a diabatic time-evolution, i.e. 
how much of the initial population remains in the first 
(diabatic) state ? 

In the mean-field approximation, the time evolution is 
given by the Gross-Pitaevskii equation 



celebrated Landau-Zener formula 
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with the nonlinear Hamiltonian 
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and e — at. The state vector is normalized to unity, 
thus the effective nonlinearity \s g = gN , where iV is 
number of particles in the condensate and g is the bare 
two particle interaction constant. Throughout this paper 
we use scaled units such that h = 1. 

The Landau-Zener transition probability is defined as 
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The original linear problem can be solved analytically 
with different approaches 0, IE IE 0- This yields the 
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for the probability of a diabatic time evolution. In the 
nonlinear case g < 0, things get quite complicated and 
the Landau-Zener probability is seriously altered. New 
nonlinear eigenstates emerge if the nonlinearity exceeds 
a critical value \g\ > gc — 2v. A loop develops at the top 
of the lowest level ^(e), while the total energy 
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shows a swallow's tail structure (cf. the left-hand side of 
Fig. The system can evolve adiabatically along this 
level only up to the end of the loop - adiabaticity breaks 
down. Gonsequently, the Landau-Zener probability does 
not vanish even in the adiabatic limit a ^ lis ,;2i| . For 
repulsive nonlinearities, g > , the situation is just the 
other way round: The loop appears in the upper level, 
thus no adiabatic evolution is possible in the upper level. 
In this paper we consider only the lower level and thus 
the attractive case g < 0- These considerations have let 
to a reformulation of the adiabatic theorem for nonlinear 
systems^based on the adiabatic theorem of classical me- 
chanics Q ■ Note also that the emergence of looped levels 
was previously studied for the quantum dimer [g. 

Several approaches were made to derive a nonlinear 
Landau-Zener formula for this problem using methods 
from classical Hamiltonian mechanics jE llol | . For subcrit- 
ical values of the nonlinearity \g\ < gc, standard meth- 
ods of classical nonadiabatic corrections yield good re- 
sults for the near-adiabatic case {ajv^ ^ 1). For the 
case of a rapid passage [v^ ja <C 1) one finds a quan- 
titative good approximation using classical perturbation 
theory with v being the small parameter for the subcrit- 
ical regime as well as for strong nonlinearities, as long as 
g < 0. Furthermore for strong nonlinearities there is a 
simple formula which provides a good approximation for 
the tunneling probability for an intermediate range of the 
parameter a. This approximation fails in the rapid limit 
as well as in the near adiabatic one. However, there is no 
valid approximation in the critical regime \g\ > gc = "i-v 
for a — *■ so far. Since one is interested in the quasiadi- 
abatic dynamic in most applications this is an important 
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FIG. 1: Total energy Q in the mean-field theory (left) and 
eigenenergies the many-particle Hamiltonian ©, (right) for 
V = 0.2, g = — 1 and A'' = 20 particles. 

deficit. Here we present a different approach wliicli yields 
good results especially in this region. 

Going back to the roots of the problem, we consider 
the original many-particle problem of an interacting two- 
mode boson field instead of the mean-field theory. We 
consider the many-particle Hamiltonian of Bose-Hubbard 
type, 

H{t) = e{t){ni - h2) + v{a\a2 + alai) 

+ f(^ii(fii-l)+fi2(r^2-l)), (6) 

where aj and aj are the bosonic annihilation and creation 

operators in the jth well and fij = ajaj is the occupation 
number operator. The eigenvalues of the Hamiltonian © 
are shown in Fig. ^ on the right-hand side in dependence 
of e for g — —1, v = 0.2 and = 20 particles. One 
recognizes the similarity to the mean-field results shown 
on the left-hand side. A series of avoided crossings with 
very small level distances is observed where the mean- 
field energy levels form the swallow's tail structure. 

For t — oo, one has e = at — > — oo and the first 
term dominates the Hamiltonian. The ground state is 
l^o) = (7V!)-i/2(aj;)^|0), where N is the fixed number 
of particles. In the spirit of the Landau-Zener problem 
we take this as the initial state for t — > — oo and consider 
the question, how many particles remain in the first well 
for t —^ +00, i.e. the effective Landau-Zener transition 
probability for the population, which is given by 

The superscripts mp and mf are introduced to distin- 
guish between the many-particle and the mean-field sys- 
tem. It will be shown that this many-particle Landau- 
Zener probability agrees well with the mean-field Landau- 
Zener probability Q. Furthermore this "back-to-the- 
roots'-procedure reduces the problem to a linear multi- 
level Landau-Zener scenario, which can be solved approx- 
imately in an independent crossing approximation. In 
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FIG. 2: The S-matrix elements |Sf,jvp in the independent 
crossing approximation (ICA) for N = 3 particles. 

this way we derive a Landau-Zener formula for an inter- 
acting BEG, which agrees well with numerical results es- 
pecially in the strongly interacting regime \g\ > gc = 2v. 

II. THE MANY-PARTICLE LANDAU-ZENER 
PROBLEM AND THE ICA 

We now consider the many-particle Landau-Zener sce- 
nario (O in detail, where the number A^ of particles is 
fixed. We expand the Hamiltonian H in the number- 
state basis \k) = [k\{N - k)\]-^/^ {al)^-''\0). Then 
the Hamiltonian is given by the matrix {£\H\k) = Hg^k 
for £, k — 0, . . . , N with the elements 

He^k = heXt) Si^k + ve {Se,k-i + Si-im) (8) 

and 

he{t) = e{t){2e -N) + ^{2f - 21N + - N) 

and the couplings V£ — {I -|- 1)(A^ — €) on the sub- 
and superdiagonal. In the Landau-Zener scenario, all 
diabatic (i.e. uncoupled) levels hgit) vary linearly in time 
as t(t) — at, however with a different offset and slope 
a{2i-N). 

As stated above, we assume that initially all particles 
are in the first well, \ip{t —oo)) = \N). Gonsequently 
one has {ni(t — oo)) = A^ and in order to derive the 
Landau-Zener probability lO we are left with the prob- 
lem to calculate {ni{t — > -|-oo)). Thus we are not inter- 
ested in the details of the time evolution. We just need 
a few elements of the S'-matrix, which is defined by 

(fc|V'(t = +oo))=^5fc£(£|V'(t = -oo)). (9) 
e 

With this definition and {£\ni\k) = kSg^k, the Landau- 
Zener transition probability Q is reduced to 

1 ^ 

^lT = ^E^I^^-^I'' (10) 

A:=0 
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FIG. 3: Landau-Zener tunneling probability in dependence 
of the parameter velocity a for v = 0.2, A'^ = 100 particles and 
different values of the interaction constant g. Numerical data 
(mean-field + and many-particle theory o) is compared with 
the ICA p2|l .fdashed line) and the resulting ICA-Landau- 
Zener formulae 1251 resp. (solid line). 



SO that only the squared modulus of the 5-matrix ele- 
ments IS'fc.ivl^ a-re of importance. 

The S-matrix elements |S'Ar_fcp are now evaluated in 
a modified independent crossings approximation (ICA, 
see appendix for details). One assumes that the system 
undergoes a series of single, independent transitions be- 
tween just two levels. The probabilities of a diabatic 
resp. adiabatic transition at a single anti-crossing are 



given pfc, AT = eyip{-TTwf. J\bk,i\) resp. qk,e = 1 - Pk,i 
according to the Landau-Zener formula l@J. Here, Wk,i 
denotes the level spacing at the anti-crossing and bk,i is 
the difference of the slopes of the two diabatic levels. The 
relevant S-matrix elements are given by 
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with the definition q^N = 1 Pnn = 0. The calcula- 
tion of the S-matrix elements by the ICA is illustrated in 
Fig. El for the case iV = 3. 

The ICA-Landau-Zener transition probability is then 
given by 
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Note that the p^^N depend on the distance between the 
levels i and N at the anti-crossing. Thus they have to 
be evaluated at different times te^N- However, the cross- 
ing time is easily calculated by evaluating ft. at (t^, at) = 
hi{ti^Pf), where hi(t) are diabatic levels as defined above. 
This yields 
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FIG. 4: Landau-Zener transition probability -Plz^(q) in the 
adiabatic limit a — ^ for i; = 0.2, g — ~1 and different 

numbers of particles: A'' = 10 (— • — ), A'^ = 20 ( ) and 

iV = 30 (— ). 



At all crossing times i^^jv, the level spacings Wi^N are 
calculated by diagonalizing the Hamiltonian matrix ||SJ|. 
As H is tridiagonal, this can be done very efficiently. 
Furthermore, the difference of the slopes is simply given 
by be,N = 2a(iV - e). 

To test this approach we compare the ICA-Landau- 
Zener formula H12|l with the Landau-Zener probability 
calculated by numerically integrating the many-particle 
Schr5dinger equation as well as the mean-field transition 
probability The results are shown in Fig. U] in de- 
pendence of a for V = 0.2, N = 100 and three different 
values of g. One observes a good agreement between the 
Landau-Zener formula (|12l) , (dashed line) and the numer- 
ical results for large g. For small values of g the ICA (|12|) 
overestimates the transition probability. These issues will 
be further discussed in section Ivl 



III. LIMITING CASES 

The linear limit g — > is analytically solvable in both 
cases. The many-particle system (jHJ reduces to the so- 
called bow-tie model, whose S-matrix was calculated in 
The mean-field dynamics reduces to the ordinary 
two-state model of Landau, Zener, Majorana and Stiick- 
elberg 0, 0, 0, Q- Not only the transition probability 
but also the whole dynamics is known exactly in terms 
of Weber functions Q • 

In the zero-coupling limit u — > the Hamiltonians be- 
come diagonal and the evolution is fully diabatic. The 
Landau-Zener transition probability tends to one. 

Most interesting is the adiabatic limit a ^ 0. As no 
subdiagonal element of the Hamiltonian matrix ^ van- 
ishes, all eigenvalues must be distinct (see, e.g. 12]). 
They may become pathologically close, but they cannot 
be degenerate. This is in fact the case: The splitting 
of the lowest levels at the anti-crossings becomes really 
small for increasing \g\. Thus all wi^n are non-zero and 
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in the extreme adiabatic limit a — > the Landau-Zener 
probabiUties p^^jv must vanish. This seems to contradict 
the mean-field results (cf. Fig. , which predicts a non- 
zero Zener tunneling probability even in the adiabatic 
limit if \g\ > gc- 

However, the parameter regime, where the ICA pre- 
dicts a vanishing Zener tunneling probability in contrast 
to the mean-field results, decreases rapidly with an in- 
creasing number of particles N. Figure 01 shows the 
Landau-Zener probability P^^^{a) for very small a, cal- 
culated within the ICA for different N, with g = gN = 
— 1 fixed. The truly adiabatic region, where Pl^^{a) « 
is negligibly small already for these quite modest num- 
bers of N. 

The mean-field theory is valid for a BEC consisting 
of a macroscopic number of atoms. In order to com- 
pare to the mean-field results we thus have to consider 
the limit of a large number of particles, N —>■ oo with 
g = gN fixed. In this macroscopic limit, the contradic- 
tion vanishes. Furthermore, this limit will prove itself as 
extremely convenient for the evaluation of Eq. H12I) , since 
all sums can be replaced by integrals which can be solved 
explicitly (cf. section lv|l . 



IV. THE MANY-PARTICLE SPECTRUM 

The only missing step towards an explicit Landau- 
Zener formula is the evaluation of the squared levels spac- 
ings w'^ j^{tk,N)- Thus one has to understand the spec- 
trum of the Hamiltonian (jHJ . We start with a discussion 
of the spectrum for e = 0, which provides an insight into 
the qualitative features which will guide us in the follow- 
ing. To keep the calculations simple, we introduce the 
operators 
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which form an angular momentum algebra with quantum 
number j = A^/2 [lUllll^. The Hamiltonian © then 
can be rewritten as 



H = 2v.L 



N 



(15) 



up to a constant term. 

In the subcritical case \g\ < 2\v\, the interaction terms 
can be treated as a small perturbation. The unper- 
turbed eigenstates are the Jz-eigenstates |j, m^) with 
niz — —J, + , j. In second order this yields the 
levels 
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FIG. 5; Spectrum of the many-particle Hamiltonian Q for 
e = 0, V — 0.2, A' = 50 particles and g = —0.1 and g = —2, 
respectively. 



up to a constant. This spectrum is illustrated in Fig. [S] 
for g — —0.1, V = 0.2 and A^ = 50 particles. The eigenen- 
ergies are nearly equidistant, with a slight increase of the 
level spacing for higher energies. 

For \g\ > 2\v\ and low energies, the interaction term 
gJ^/N^ dominates the Hamiltonian. The eigenstates 
with quantum numbers \j,±mx) are doubly degenerate 
with eigenenergy — gm'^/N'^. The perturbation 

2vJz removes this degeneracy only in the 2|ma;|-th or- 
der. Thus the low energy eigenstates (corresponding to 
the high \mx\ states) appear in nearly degenerate pairs. 
However this approach fails if the energy scale of the per- 
turbation 2vJz becomes comparable to the unperturbed 
eigenenergy. Estimating the energy scale of the pertur- 
bation as £^niax/2 = 2vj/2, perturbation theory fails for 
\g\m1/N^ < vN/2. Instead, Bogohubov theory provides 
the appropriate description for the high energy part of 
the spectrum. We are dealing with an attractive inter- 
action 5 < 0, so that the highest state in the mean-field 
approximation is the state with equal population in the 
two modes. So the standard Bogoliubov approach is valid 
for the highest state instead of the ground state. One 
finds that the high energy part of the spectrum is given 
by_£'„ — En — lu(N — n) with the Bogoliubov frequency 
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To clarify this issue, the spectrum is plotted in Fig.|31for 
g ~ —2, V — 0.2 and A^ = 50 particles. One clearly sees 
the nearly degenerate pairs of eigenvalues for low energies 
and the approximately equal spacing of the high-energy 
eigenvalues. The distance of the two highest levels is 
given by the Bogoliubov frequency (|17l) . 

Now we come back to the squared level splittings 
ui^ Ar(ifc,Af): beginning with the supercritical regime \g\ > 
2v. Figure |5| shows an example of the squared level 
splitting for v = 0.2, TV — 100 particles and g = —0.1 
resp. g = —1. Later, we consider the macroscopic limit 
N — > 00, g — > with g — gN fixed. For this issue we 
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FIG. 6: Squared level splitting w^{x) in dependence of the 
scaled index x = l/N for v — 0.2 and g — —0.1 resp. g = 
— 1. Numerical results (o) are compared to the approximate 
formulae I19II resp. CT . (solid lines). 



plot the squared level splittings versus the rescaled index 
X :— £/N e [0, 1]. With increasing N, the curve plotted 
in Fig. 1^1 remains the same, only the actual points move 
closer together. Thus one obtains a continuous function 
w'^{x) in the limit N oo. 

As argued above for e = 0, the lower levels appear in 
approximately degenerate pairs. By the same arguments 
one concludes that this is also true for for the first level 



crossings. Thus, 



is effectively zero for £ < £c resp. 



X < Xc- The critical index Xc can be estimated as de- 
scribed above for e = 0. It is found that this estimate 
gives the correct results up to a numerical factor a of 
order 1. Thus we conclude that 



« 1 — a^2v/\g\. 



(18) 



A very good agreement of this formula to the numerical 
results was found for a = 1.14. 

For X > Xc the squared splittings increase approxi- 
mately linear. In the high energy limit corresponding 
to X ^ 1, the level splitting is given by the Bogoliubov 
frequency introduced above. In conclusion, the squared 
level spacing can be approximated by 



w'^{x) 



to 
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'-H{x - Xc), 



(19) 



where H{x — Xc) denotes Heaviside's step function. 

In the subcritical regime \g\ < 2v, one can use the 
results from perturbation theory described above (cf. 
Eq. lltill ). At time tg^^ — ~g£/2a one must evaluate 
the level splitting S^^j+i — Eg^j (note that the levels are 
labeled by — ~j, ~j + 1, . . . ,j with j = N/2). Again 
we consider the limit — > oo with g = gN fixed. After 
a little algebra one finds that the relevant level splitting 
is in linear order given by 
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16w2 + 4g-3gV4 
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(20) 



in terms of the scaled index x = £/N . 

The approximate results for the squared level splitting 
w'^{x) for \g\ > gc, Eq. lO, and for \g\ < gc, Eq. ^U^, 
are compared with the numerical results for A^ = 100 
particles in Fig. El One observes a good agreement. 



V. TOWARDS AN EXPLICIT LANDAU-ZENER 
FORMULA 

Using the formulae for the squared level splitting de- 
rived in the previous section, the ICA-Landau-Zener 
transition probability (|12|l can now be evaluated explic- 
itly. In the spirit of the of the macroscopic limit N ^ oo, 
the sums are replaced by integrals according to 
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The difference of the slopes bi^N — 2a{N — £), which en- 
ters the formula is also rewritten in terms of the rescaled 
index x = £/N: 



bLN -> 2aA^(l - x) =: Nb{x) 
Thus one finds 
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In the supercritical regime \g\ > gc we start by evalu- 
ating the integral over x in Eq. (|23|) . Substituting w'^{x) 
and b{x) from Eq. H19(l and H22() and carrying out the 
integral yields 
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for y > Xc and zero otherwise. The Landau-Zener tran- 
sition probability (|23|l is then given by 
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(25) 



with the abbreviation u = 7rw^/2a and Xc defined 
in Eq. 11811. Here, 7 denotes the incomplete gamma- 
function [131 . 

In the subcritical regime \g\ < gc, one finds by substi- 
tuting Eq. (|23|l into (|23|) . that the Landau-Zener transi- 
tion probability is given by 
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FIG. 7: Landau-Zener tunneling probability in dependence of 
the interaction constant g for a parameter velocity a = 0.01. 
Numerical data (mean-field + and many-particle theory o) is 
compared to the ICA-Landau-Zener formula 1251 1. fsolid line) 
for V = 0.2. 



with the abbreviations cq = 7r(wo + 2wqWi) /2a and ci = 
■kwqWi/ a. To keep the calculations feasible, we kept only 
terms linear in x resp. y in the exponent consistent with 
Eq. 

To test the validity of our approach we compare the 
ICA-Landau-Zener formulae H25I) and (|2(j|l to numeri- 
cal results obtained by integrating the SchrSdinger equa- 
tion for mean-field Hamiltonian (0) as well as the many- 
particle Hamiltonian 10. The Landau-Zener tunneling 
probability in dependence of the interaction constant g 
is plotted in Fig. [T] for a = 0.01 in dependence of the 
velocity parameter a for different values of g in Fig. 

One observes a good agreement of the ICA-Landau- 
Zener formula with the numerical results in the critical 
regime \g\ > Qc- Especially the increase of the tunnel- 
ing probability with increasing \g\ for small a is well de- 
scribed by our model. This problem could not be solved 
with previous approaches 9,10]. The approximation gets 
worse for larger values of a because the assumption that 
the Zener transitions are well separated becomes doubt- 
ful for such a large parameter velocity. The ICA thus 
underestimates the tunneling probability. 

In the subcritical case \g\ < gc, the proposed ICA- 
Landau-Zener formula does not work as well. In fact 
the tunneling probability is overestimated for small a be- 
cause the ICA itself is not a very good approximation in 
this case. The adiabatic levels do not show well sepa- 
rated avoided crossings, instead the levels splittings are 
nearly constant over a long interval of the parameter e. 
For larger values of a one faces the same problems as 
in the supercritical case and the tunneling probability is 
underestimated. Another ansatz, using e.g. perturbation 
theory with respect to the solution of the noninteract- 
ing problem should be better suited to this problem. 
Note however, that the deviations are mainly due to the 
ICA itself and to the approximation of w'^{x) made in 
this section. 



VI. CONCLUSION AND OUTLOOK 

In conclusion, we have derived a Landau-Zener for- 
mula for an interacting Bose-Einstein condensate from 
first principles. To this end we considered the original 
two-mode many-particle Landau-Zener scenario. It was 
shown that the resulting Landau-Zener formula agrees 
well with the numerical results calculated for the many- 
particle problem as well as within the mean-field approx- 
imation. 

In the future, it would be of interest to relate our cal- 
culations to the respective problem in the Heisenberg pic- 
tures. Here, complex eigenfrequencies may occur for the 
dynamics of the creation/annihilation operators, leading 
to spontaneous production of quasi-particles and hence a 
dynamical instability. For the non-interacting case, this 
problem has been solved analytically 18]. 

Another issue is the discussion of nonlinear Landau- 
Zener problems for more than two levels. First results 
for three level system were reported only recently [l9|. 



Appendix: The independent crossings 
approximation 

Let us first briefly recall the dynamics of a two-level 
Landau-Zener system described by the Hamiltonian 

The diabatic and adiabatic energy curves are plotted in 
Fig. m The S-Matrix in the sense of Eq. |^is given by 

S^(^^l) (2S) 

with p — exp (— TTw^/]/?! — (32\) and q ~ — p^ . The 
probability of a diabatic passage is therefore given by 

the Landau-Zener-formula Flz = p^ = exp ^— j^^^^^, 

and the adiabatic transition probability by 1 — Plz — <^ ■ 
They depend only on the relative slope of the diabatic 
levels |/3i — /32l and the coupling which is equivalent 
to half of the gap between the adiabatic energy-levels at 
the avoided crossing. 

The simplicity of the solution of the two-level system 
and the observation that the transitions between two adi- 
abatic levels in a multilevel Landau-Zener system takes 
place only in a very narrow region around the crossing of 
the two corresponding diabatic levels leads to a simple ap- 
proximation. If all crossings are well separated they can 
be considered as independent of each other and each of 
them is described by the two-level-Landau-Zener model 
where the couplings between the relevant diabatic lev- 
els are given by the nondiagonal terms of the Hamilto- 
nian. This approach is called the "independent crossing 
approximation" (ICA) in the literature. It is of great 
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FIG. 8: Diabatic (dash-dotted line) and adiabatic (solid line) 
energy levels of the two-level Landau-Zener model 1271 . 

importance for the study of multilevel Landau-Zener dy- 
namics because of a surprising feature: The ICA turns 
out to give the exact results for all known exactly solvable 
multi- level Landau-Zener scenarios ,11.,;2Q3. Furthermore 
it has been shown that the ICA always gives the correct 
results for the diagonal S'-matrix elements with minimal 
and maximal slope 0, 0| . 

Of course there are also examples where the ICA fails, 
as for example for the simple three level Hamiltonian 

(at ^ a V w \ 
V Q ■ (29) 
w —at + a I 

The adiabatic and diabatic levels are plotted in figure 
El The diabatic transition probability for the third level 
is exactly given by the ICA. But if we look at the 
S-Matrix element ^32 we find that the ICA predicts it to 
be zero, because the coupling matrix element vanishes, 
(2|iJ|3) = 0, independent of a and a, which isn't true. 
The second and the third diabatic levels do not couple 
directly, but for not too large values of a the indirect 
coupling via the second diabatic level can't be neglected. 
This coupling manifests itself as an avoided crossing be- 
tween the two adiabatic levels, which turns into a real 
crossing only in the limits a 00 and a ^ 0. For finite 
values of a the transition probability between the third 
and the second diabatic levels is small but nonzero. 



To get a better approximation one should recall the 
two level system, where the coupling between two dia- 
batic levels is equivalent to half of the level-splitting of 
the corresponding adiabatic levels. Therefore one can use 
a modified ICA where the couplings are not given by the 
nondiagonal elements of the Hamiltonian but half of the 
level splitting between the relevant adiabatic levels. This 
approximation doesn't inherit the benefit of providing 
the exact results in the special cases where the original 
ICA did, but provides a good approximation even in the 
cases where the ICA fails. Therefore it is better suited 
for our purposes. The performance of the approximation 
is limited by the fact that the single avoided crossings 
must be well separated so that the transition regimes do 
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FIG. 9: Diabatic (dash-dotted line) and adiabatic (solid line) 
energy levels of the three-level Landau-Zener model 1291 for 
a = 0.2, a = 0.5, v = 0.2 and w = 0.3. 



not overlap. In the present case this is improved with in- 
creasing nonlinearity. Note that to simplify matters this 
modified ICA is denoted as ICA throughout the paper. 
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